Solitary waves in mixtures of Bose gases confined in annular traps 
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A two-component Bose-Einstein condensate that is confined in a one-dimensional ring potential 
supports solitary-wave solutions, which we evaluate analytically. The derived solutions are shown 
to be unique. The corresponding dispersion relation that generalizes the case of a single-component 
system shows interesting features. 

PACS numbers: 05.30.Jp, 03.75.Lm 



Introduction. Cold atoms provide an ideal system for 
the study of nonlinear effects including solitary waves. 
The study of solitary waves in these systems has a num- 
ber of interesting features that extend the results of the 
well-known nonlinear Schrodinger equation in homoge- 
neous systems, first studied by Zakharov and Shabat [l| 
three decades ago. 

Among the many novel aspects associated with the 
physics of solitary waves in trapped gases 0,0, IH we note 
(i) the presence of an external trapping potential which 
renders these systems finite, (ii) the fact that they can 
be constructed as quasi-one or quasi-two dimensional sys- 
tems, and (iii) the possibility of bound states of solitary 
waves in two- (or multi-) component gases. We refer the 
reader to the article by Carretero-Gonzalez et al. @ for 
a review of the extensive work that has been performed 
on all these problems. 

Remarkably, recent experimental advances now per- 
mit the realization of elongated, quasi-one-dimensional 
traps. Solitary waves have been observed in elongated 
traps as well as in more spherical traps 0|. More 
recently, solitary waves have also been observed Q in a 
two-component Bose-Einstein condensate. 

In the present study we examine the problem of soli- 
tary waves in a two-component Bose-Einstein condensate 
B 0, II H M, El El El H on an infinite line. Making 
an ansatz for the density of the two species, we derive 
exact solutions of the two coupled nonlinear equations, 
which describe the order parameters of the two gases and 
derive the profiles of grey-grey and grey-bright solitary 
waves, with the assumption that the density is constant, 
either zero or nonzero, far from the center of the two 
waves. We then demonstrate the uniqueness of the de- 
rived solutions. Finally, we derive the dispersion relation 
associated with these solitary waves, imposing periodic 
boundary conditions, i.e., assuming that the bosons are 
confined on a ring. Since the derived solutions are ex- 
ponentially localized, the solutions with periodic bound- 
ary conditions can be approximated accurately by the 
present infinite-line calculations if the ring is sufficiently 
large. We compare the resulting dispersion relation with 
that of a single-component Bose gas confined in one di- 
mension under periodic boundary conditions and inter- 
acting via a contact potential, that was first studied by 



Elliott Lieb [151 . 

Model. A Bosc-Einstcin condensate consisting of two 
distinguishable species, A and B, is described within 
the mean-field approximation by the two order param- 
eters ^a(x, t) and &b(x, t), which satisfy two Gross- 
Pitaevskii-like equations, i.e., two coupled nonlinear 
Schrodinger equations. Assuming for simplicity equal 
masses for the two species, Ma = Mb = M, and set- 
ting h = 2M = 1, these equations have the form 



id^j/dt = -dHj/dx 2 + 7(|*j| 2 + \^ k \ 2 )^j 



(1) 



where we have also assumed equal s-wave scattering 
lengths oaa = clbb = &ab = a for all elastic collisions 
between the atoms. The coefficient 7 that multiplies the 
nonlinear terms is thus equal for all terms and is also 
oc a. In the above equations we have also chosen the nor- 
malization condition j \^j\ 2 dx — Nj , where Nj is the 
number of atoms in each component. 

Solitary waves. We first determine the solitary wave 
profiles that result from Eqs. ([1]). Writing ^j(x, t) = 

rij (x, t) e 1 *' ( x '*) and separating real and imaginary 
parts, we get two continuity equations for the two species, 
as well as two Euler equations, namely 



drij/dt 



d^/dt = (VnJ) 



2(n^.)' 



0, (2) 
j(nj+n k ), (3) 



where the primes denote spatial derivatives. From the 
asymptotic behavior of Eqs. |T]), we find that cx e~ vt , 
where /1 = ~f(n A + n B ) is the chemical potential and n° 
is the density of species j at \z\ = ±00. 

To derive solitonic solutions, we assume a constant and 
common velocity of propagation u for the two waves, and 
therefore we assume that &j(x,t) = \J n~j (z)e l ^ j ( z )~ ^'1 , 
where z — x — ut. We note that while the density 
rij = \^j(x, i)| 2 is indeed a function of z = x — ut only, 
the phase has a more general dependence on x and t as 
indicated. 

Grey-grey solitary waves. The assumption of 
travelling-wave solutions allows us to convert the time 
derivatives into spatial derivatives. The continuity equa- 
tions may then be integrated to give 



0' =(«/2)(l-n> J -)- 



(4) 
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Since the local fluid velocity is Vj = Tuft^/M, we see that 
u(l — nS/rij). The Euler equations now take the 



form 



u 

T 



(n ) 2 \ 

- 1 + 7(»J ~ »° + »* - (5) 



where we have imposed the boundary conditions nj — > n° 

and (y/rij) — > for |z| — > oo. 

To integrate the above equations, we make the follow- 
ing ansatz, which is also consistent with the boundary 
conditions, 



n B - n° B = n(n A - n° A ). 



(6) 



This ansatz allows us to perform the integration and ob- 
tain 



nl/2 = + K )n A - u 2 /2}(n A - n° A ) 



^2 



n'i/2 



h- 



\2 



-n B - U /2](n B - n B ) 



(7) 
(8) 



The consistency of Eqs. ©, 0, and (jHJ) requires that 
tib = ^n^, which establishes the value of re as re = n B /n A 
independent of the velocity of propagation u. Since the 
sign of re is clearly positive, there can be no dark-antidark 
solitonic solutions that are consistent with the ansatz of 
Eq. ([6]). In other words, if the are nonzero, the waves 
must both be density depressions, i.e., both must be grey 
solitary waves. 

The solution of the above equations in the case of grey- 
grey solitary waves has the form 



-rij cos 2 



6 /cosh 2 (z cos 6/£i 



(9) 



where sin 2 9 — u 2 /c 2 with c 2 = 277^ (1 + re) is the (com- 
mon) speed of sound and where the coherence length 
is given by l/^ 2 , = 70,^(1 + re)/2. 

Grey-bright solitary waves. The above analysis re- 
quires some modifications in the case where the density 
of one of the components vanishes at infinity. Assuming, 
for example, that n A — ► n A ^ 0, and jib — > 0, the phase 
of A is still given by Eq. (U|). However, the local velocity 
of species B is constant and equal to it, since <j>' B = u/2. 

The equations obtained in this case are 



u 

T 



(^T-- 1 ) +~f(n A ~n A + n B ), (10) 



A 

(V^b) 



= l{n A 



n B ) 



(11) 



where s is the limit of (■y/n B )" / \JnB as \z\ — > 00. The 
parameter n B in the case of grey-grey solutions is now 
replaced by s, which can be determined from the number 
of particles Nb in the bright component B. The ansatz 
for the solution assumes the form tib = n(n A — n A ) with 
re negative. Here, re depends on the propagation velocity 
U, in contrast to the result obtained above for grey-grey 
solitons. 



The integration of Eqs. (fP0"|) and (fTTj) yields 



\2 



,'2 



f/2 = + K )n A - u 2 /2}(n A - n A ) 

K -f- 1 

12 = [7 n B + 2s}n 2 B . 



(12) 
(13) 



Consistency of these equations with the ansatz demands 
that 

M 2 /2 = 7 (re+l)n^-2s, (14) 

and therefore re > — 1. Since re must also be negative, 
-1 < re < 0. 

The solution of Eq. (fl2|) for the grey component is the 
same as in Eq. © . The solution of Eq. (|13[) for the bright 
component is 



1 



11 B 



— Ml a COS 



2£b cosh {z/£b) cosh (z cos6>/££>) 



(15) 



where J ns dz = Nb, in accordance with the normaliza- 
tion condition, and s = l/£fj = cos 2 0/£,fj. 

For completeness we mention that on a ring of finite 
radius, our equations also support bright solitary waves 
in both components provided that 7 is negative. We will 
examine this problem in a future study. 

To get some insight into the above results, it is instruc- 
tive to write the initial equations taking into account the 
ansatz, 



.d^ A 



d 2 ^ A 
dx 2 
d 2 ^ 



+ 7[tu(1 + re) + n% - Kn° A ]V A , 



(16) 



+ l[nB- 



1 



11 B 

re 



*s. (17) 



dt dx 2 re 

In a sense, the ansatz decouples the two equations, al- 
though consistency forces additional conditions as shown 
above. The final terms on the right of Eqs. (fTS)) and (fT7|) 
either vanish for the grey-grey case, or are constant for 
the grey-bright case and are not important. The impor- 
tant terms are those involving the combinations n A (l + K.) 
and 71b(1+1/k) or, equivalently, the "effective" couplings 
U A = 7(1 + re) and U B = 7(1 + 1/re). Clearly U A = kU B - 
In the case of grey-grey solitary waves, U A , Ub and re 
are all positive, which is consistent with the fact that 
the solitary waves are grey in both components. On the 
other hand, in the case of grey-bright solitary waves, Ub 
is negative and U A is positive since — 1 < re < 0. This is 
consistent with a bright wave in the B component and a 
grey wave in the A component. 

Uniqueness of the solutions. It is useful to regard the 
present problem as that of the motion of a particle in a 
two-dimensional potential, with ^JfTj playing the role of 
spatial coordinates and z that of time. For example, in 
the case of two grey solitary waves, 



{d^/dz) 2 /2 
where 



/^E/dz) 2 /2 + V(n A , n B ) = 0,(18) 



V 



(n A 



,0 \2 



(nB 



11 



n A 



n° A 



n B 

n B - 



(19) 
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In the case of grey-bright solutions the potential becomes, 



V = 



u 2 (n A - n° A ) 2 
8 n A 



7 



1 



{n A -n° A + n B ) 2 - -sn B - (20) 



In both expressions for V above, we have imposed the 
boundary condition that V vanishes when n A and Ub 
have their asymptotic values. 

In the grey-grey case, the directional derivative of the 
potential perpendicular to the straight line defined by the 
ansatz vanishes, 



(a 



dV dV 



dJn A ' dJriB 



= 0. 



(21) 



Starting at the asymptotic field values at "time" z = 
— do, the system will move along this minimum in V and 
return to its starting point at z = +oo. 

The solutions found above for the specific boundary 
conditions are unique. The only physically interesting 
trajectory must start at rest on the V = contour with 



A a ^ V "■B 

same point. Linearization of the two coupled nonlinear 
equations for \z\ — > oo sets the asymptotic form of the two 
solutions. Given this asymptotic form, the solutions fol- 
low a well-defined path for all z determined by the poten- 
tial of Eqs. (fl9|) and (f20|) . In other words, there is one and 
only one trajectory which will result from starting the 



system at the point with (y/nX, y/n B ) = ( 



a- 



Dispersion relation. Having calculated the profiles of 
the solitary waves, it is instructive to evaluate the cor- 
responding dispersion relation. In the case of a single- 
component Bose gas in one spatial dimension, Eliott 
Lieb [151 ] found that when bosons interact via a con- 
tact potential, the excitation spectrum consists of two 
branches. The one corresponds to the usual Bogoliubov 
mode, while the other was later shown by Kulish et al. 
and by Ishikawa and Takayama to correspond to solitary 
waves [HG1 EE El- 

The present calculation generalizes these results to the 
case of a two-component system with some similarities 
and some differences. Starting with grey-grey solitary 
waves, the kinetic energy is given by 



KE 4 



[(dy/nj/dz) 2 + midfa/dz) 2 ] dz. (22) 



From Eqs. (j4|) and ([8]) the total kinetic energy is 

KE = KE A + KE B = (1 + k) 2 | J {n A - n° A ) 2 dz.{23) 

The interaction (free) energy less the infinite energy of 
the background density of the two components is 



PE - ji A N A - ti B N B = | / (ji A -n° A + n B - n° B ) 2 dz, 



(24) 



which is equal to KE. The total energy is thus 

E/E GG = (4/3)(l + k) 1 ' 2 cos 3 0, (25) 
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FIG. 1: The dispersion relation E = E(P) (solid line), dE/dP 
(dashed line), and u (dashed line) as functions of P for the 
case of grey-grey (upper graph) and grey-bright (lower graph) 
solitary waves. The energy is measured in units of Egg in the 
grey-grey case, and in Egb in the bright case. The derivative 
dE/dP and u are measured in units ^2n A ^y. The figures show 
the case n A = n° B in the grey-grey case, and Nb is set equal 
to n A £o in the grey-bright case. The equality dE/dP — u is 
satisfied in both cases. 



where E GG = jn° A (l + K)(n A £ ) with l/£g = i2° A -//2, or 
E GG = (1 + «)[27(n^) 3 ] 1 / 2 . Turning to the momentum, 



P/n° A = (1 +k)(tt- 20 -sin 20), 



(26) 



where we have imposed periodic boundary conditions; 
for a detailed analysis of the derivation of Eq. (f2"6")h see 
Ref. [H[. When Eq. (J26J) is combined with the result of 
Eq. (|25|) . we get the dispersion relation E = E(P) plotted 
in the upper graph of Fig. 1 for the case n A = n° B . 

The dispersion relation for grey-bright solitary waves 
is calculated in a similar way. The kinetic energy of the 
two components is 



KE A = ( 7 /2)(k + 1) l(n A -n° A ) 2 dz 



(27) 



KE B = i«7)AT B (l + K) + |«(K+l) J(n A -n° A ) 2 dz, 

(28) 

while the potential energy minus the infinite energy due 
to the background density of the A component is 

PE-fi A N A = ^{l + K ) 2 [(n A -n° A ) 2 dz. (29) 
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Collecting terms yields 

E/E GB = 1(1 + k) (l-l^cos 2 ^ , (30) 

where Eqb — N B n° A j_. The momentum of the two com- 
ponents is given as 

P A /n° A = tt- 26» -sin 29, (31) 
P B = N B sin 0y/T+Hy/n A >y/2. (32) 

The parameter k is now a function of u (and therefore of 
9) . Their relationship is established by the normalization 
condition, 

N B /(n A to) = -2k cos 6>/\/TT^. (33) 

From Eqs. J3T), (32]) and {33]) we find that P = P 4 + P B 
is given by 

P/n° A = tt - 20 - (1 + k) sin 29. (34) 

Combining Eqs. (|3"D|) , (f3"3"|) and (|3"4"]) . we obtain the disper- 
sion relation for a grey-bright solitary wave, which is plot- 
ted in the lower graph of Fig. 1 for the case n A £o = N B . 
Remarkably, the equation u = dE/dP is satisfied in both 
graphs. This justifies the picture of a solitary wave - or 
rather a bound state of two solitary waves in the specific 
case - as a particle. 

While the dispersion relation shown in the upper panel 
of Fig. 1 is qualitatively the same as in the case of a 
single-component system, with a constant slope for long- 
wavelength excitation (corresponding to sound waves) 
and zero slope at the maximum value of the momentum 
(corresponding to dark solitary waves in both compo- 
nents), the lower panel is different for long- wavelength 
excitations. In this limit, the dispersion relation has a 
vanishing slope. The limit of a static dark solitary wave 



is qualitatively the same in the two cases; small ampli- 
tude sound waves show a different behaviour. 

Specifically, for the case of grey-grey solitary waves 
when P — > 0, we find E = cP, with the speed of sound 
c = ^/2~/(n A + n B ) identical to the value obtained by 
linearizing the original equations [2(|. For the case of 
grey-bright solitary waves when P — > 0, the dispersion 
relation vanishes quadratically, E — > P 2 /N B . Then, 
dE/dP — > ^/2n^7(l + k), which agrees with the value 
of u obtained from Eq. (fl"4"|) for s — ► 0. 

Conclusions. Solitary waves in a two-component Bose- 
Einstein condensate confined in a ring potential show 
interesting physics. Considering repulsive effective in- 
teractions, this system supports bound states of grey- 
grey and grey-bright solitary-wave solutions, which are 
unique. The dispersion relation resembles that found for 
single-component solitary waves for large density depres- 
sions for both grey-grey and grey-bright solitary waves. 
For small density depressions, i.e., for sound waves, there 
is a significant difference between grey-grey and grey- 
bright waves. For long wavelength excitation the velocity 
of propagation u is maximum in the first case, which is 
equal to the speed of sound, and it decreases monotoni- 
cally down to zero when the grey wave become dark, as 
in the case of a single component. On the contrary, in 
the case of grey-bright waves, u vanishes for small den- 
sity variations, i.e., for long- wavelength excitations, and 
it vanishes also when the grey component becomes dark, 
having a maximum value in between, as seen in Fig. 1. 

Experimental verification of the derived dispersion re- 
lation shown in Fig. 1 would confirm the expected spec- 
trum for the first time, as this has never been confirmed 
in any nonlinear system so far. Such an experiment 
should be possible with use of the method of Bragg spec- 
troscopy [HI, H2, Hil > which has been developed to probe 
the excitation spectrum of cold atomic systems. 
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